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We describe a self consistent magnetic tight-binding theory based in an expansion of the 
Hohenberg-Kohn density functional to second order, about a non spin polarised reference density. 
We show how a first order expansion about a density having a trial input magnetic moment leads 
to the Stoner-Slater rigid band model. We employ a simple set of tight-binding parameters that 
accurately describes electronic structure and energetics, and show these to be transferable between 
first row transition metals and their alloys. We make a number of calculations of the electronic 
structure of dilute Cr impurities in Fe which we compare with results using the local spin density 
approximation. The rigid band model provides a powerful means for interpreting complex magnetic 
configurations in alloys; using this approach we are able to advance a simple and readily understood 
explanation for the observed anomaly in the enthalpy of mixing. 
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I. INTRODUCTION 



There is much subtlety connected with itinerant mag- 
netism in transition metals that one would neverthe- 
less wish to capture in a simple model. Recently 
an interatomic potential including magnetism has been 
proposed^ which will prove very useful for molecular dy- 
namics, but will not be able to describe electronic struc- 
ture effects such as the competition between ferro- and 
antiferromagnetism, or the sudden collapse of the mo- 
ment in hcp-Fe under pressurc^^i^ There are very much 
greater difficulties attendant on interatomic potentials 
employing a term in the energy which is linear in the 
magnetic moment. Almost certainly a minimum re- 
quirement of a simple model is that it contains an ex- 
plicit account of the electron kinetic energy. This is be- 
cause inter-site magnetic interactions are carried by the 
hopping matrix elements of the one-electron part of the 
Hamiltonian, not by inter-site two-electron Coulomb in- 
tegrals, and so Heisenberg and Ising models are not ap- 
propriate to discuss itinerant magnetismp^i^ The tight- 
binding approximation on the other hand provides just 
such a descriptiouj^iiSiiii in its most economical form it 
becomes a bond order potential recently described for 
transition metals by Drautz and Pettifor.'* Whether a 
magnetic bond order potential will appear remains to 
be seen; as we find below and as pointed out in Q 
an accurate prediction of some magnetic affects requires 
quite detailed structure in the density of states near the 
Fermi level. Magnetic tight-binding has been proposed 
many times using two slightly different self consistent 
schemes. The firslj^i^ii^ is based in a rigid band approx- 
imation first used by Andersen et aL— i^i^ in the con- 
text of the local spin density approximation (LSDA)!^^ 
The non spin polarised density of states is allowed to 
split rigidly as a result of on-site exchange and correla- 



tion interactions and an energy functional (equation ([7]) 
below) is minimised. This procedure may also be used 
in atomistic simulation if applied to the local density of 
states site by site; and provides a simple way to include 
effects such as magnetic pressure at crystal defects and 
site dependent magnetic momentsjiii A second more gen- 
eral approach is a self consistent scheme in which the 
rigid band approximation is lifted and both the density 
of states and the exchange splitting are determined self 
consistently^ We are motivated to recast this procedure 
into our recently proposed self consistent polarisable ion 
tight-binding model j^^i^^'^^ based on an expansion of the 
Hohenberg-Kohn functional^S to second order in a ref- 
erence electron density. We will employ a non spin po- 
larised input density, which may seem surprising but is 
consistent with the Stoner form of the LSDA which ex- 
pands the exchange correlation potential to linear order 
in the magnetic moment J^ii^ii^ 

Having described magnetic tight-binding from the point 
of view of the second order expansion, we construct very 
simple tight binding models for Cr, Fe and Co which 
we expect to be transferable to other transition met- 
als and their alloys. Finally we address an outstanding 
question in the thermodynamics of Fe-Cr alloys, namely 
the anomalous negative enthalpy of mixing at the Fe- 
rich end of the phase diagram;^ It is now well known 
that whereas over most of the concentration range Fe 
and Cr are immiscible, at low concentrations Cr is sol- 
uble in Fe, with a negative enthalpy of mixing. An ex- 
planation based on a phenomenological Ising model has 
been proposed;^ and a classical potential has been fit- 
ted to reproduce the phase diagramj^i Recent LSDA 
calculations^^ revealed that Cr atoms favour cluster- 
ing except at low concentrations when there is a repul- 
sive interaction between Cr impurities. Klaver et al^ 
pointed to this repulsive interaction in order to explain 
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the negative to positive upturn in the enthalpy of mix- 
ing at concentrations in the range 8-12 atomic percent 
Cr. Bandstructure arguments have been put forward 
based on densities of states within the coherent poten- 
tial approximation,.?- but these were rather far removed 
from the actual densities, somewhat invalidating the con- 
clusions. We are able to advance explanations for these 
phenomena using tight-binding calculations which are re- 
markably close to our LSDA results and which give rise 
to a ready explanation easily understood within the rigid 
band Stoner-Slater picture of itinerant magnetism. 

The structure of the paper is as follows. In section |TT] 
we describe how to include spin polarisation into the self 
consistent polarisable ion tight-binding model; and we 
decribe how the rigid band Stoner-Slater picture may be 
recovered from the same framework in section IIIIl In 
section IIVI we deduce parameters for a simple, transfer- 
able, non orthogonal tight-binding model for transition 
metals. We apply this model to pure Fe and Cr in sec- 
tion |V] and to Co in section IVII In section IVIII we apply 
the model to structural energetics of pure Fe. In sec- 
tion IVIIII we address the electronic structure of Fe-Cr 
alloys and in section |IX] describe the use of the self con- 
sistent rigid band model to predict the magnetic struc- 
ture and energy. We propose an explanation of the en- 
thalpy anomaly in section |Xl and conclude in section IXll 
In Appendix [X] we show how an equivalent form of the 
electron-electron interaction energy to that derived in 
section |TT] may be obtained from a multiband Hubbard 
model as used in LDA-I-J7 theory, which exposes the ne- 
glect of self interaction correction in LSDA and our mag- 
netic tight-binding while indicating how this could be 
put back into a tight-binding scheme. In Appendix |B] 
we describe non orthogonal self consistent tight-binding; 
in particular we show that in this case self consistency 
leads to adjustment of the hopping integrals in addition 
to the on-site increments, and we illustrate the origin of 
additional contributions to the interatomic force arising 
from bond charges. 



charge densities. Hq is then the Hamiltonian whose effec- 
tive potential is generated by Pi„.— We introduce a spin 
density p = J2a Tr/5°" = J2a P'^ = + P j the electron 
spin taking the value a — ±1 in units of ^h. Minimi- 
sation of the Hohenberg-Kohn functional leads to two 
Kohn-Sham equations,— in atomic Rydberg units, 

(-v2 + y-) r = er 

in an effective potential 



where V-a is the Hartree potential, V^^t 
tential due to the ions and 



VZ = 



5E^. 
If? 



the external po- 



(1) 



is the exchange and correlation potential. In the absence 
of a magnetic field (which we could include as a Zeeman 
term in V^^t) this is the only term which is spin depen- 
dent. The corresponding Hohenberg-Kohn-Sham energy 
functional is (we may supress the symbol dr under an in- 
tegral sign) 



E^+E^^ [p+,p-] +E^ 



in which T is the kinetic energy operator, E^ is the 
Hartree energy and E^z is the ion-ion interaction. This 
is expanded about the reference non spin polarised den- 
sities 



II. SELF CONSISTENT TIGHT-BINDING 
INCLUDING MAGNETISM 



and we define 



In our self consistent polarisable ion tight-binding model 
we express the electron Hamiltonian as 

H = Ho + H'. 

The first term is the usual non self consistent tight- 
binding Hamiltonian of non interacting electrons.^ H' de- 
scribes electron-electron interactions and is constructed 
so as to represent second order terms in the expansion 
of the Hohenberg-Kohn density functional about a refer- 
ence density Pi^r^ We take it that Pi„ is constructed by 
overlapping spherical, neutral, non spin polarised atomic 



Sp" ^ P" - P^; Sp = p- p,,, = Sp+ + Sp-. 

The exchange and correlation energy is expanded to sec- 
ond order in Sp" to give 

E^^ [p+,p-]=EZ + Y^ J V-Sp'^ 
2^ J J 6p^6p^' 
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The Hohenberg-Kohn total energy, exact apart from the 
neglect of terms higher than second order in E^^ is^ 



- Idr Idr' > "2 



Sp{r)Sp{r') 



Sp'^{r)6p'^'{r') 



Tj-rM (r') ■ (2) 



the total density and the magnetic moment (equation ([5]) 
below). First, we note that the magnetic moment m is 

m ^ p^ — p^ = 6p^ — (5/9^ = 6m 

since the input density is non spin polarised. We then 
find, using © and d!]) 



Sp"^ 



whereas^i2^ 



/ = -2 



S^E^. 



The first two lines amount to the Harris-Foulkes 
functional liii^Si^ The second line is represented by a 
pairwise repulsive energy, E^^^^ , in the usual tight-binding 
models. In our self consistent polarisable ion tight- 
binding model we approximate the third line as the elec- 
trostatic interaction energy between point multipole mo- 
ments of the charge transfer. The fourth line is the exten- 
sion of the on-site electron-electron interaction Hubbard 
term to the spin polarised case, and we now examine this 
term in more detail using ([T]) by writing 

<7<7' 

Here we have supressed the r-dependence, firstly be- 
cause all off-diagonal Coulomb terms are relegated to 
the Madelung energy (the third line in equation ^) 
in our tight-binding model, recognising that itinerant 
magnetism is a consequence of on-site exchange and 
correlation;^ and secondly because in our tight-binding 
model we will be using a local orbital basis to represent 
the spin density. 



The quantity 



Sp (5/9 ^ 



= u 



(3) 



is the direct Coulomb, correlation only, interaction 
strength between unlike spins described by the Hubbard 
U parameter. On the other hand the quantity 



sv- 



5p^ Sp 



U -I 



(4) 



reflects the lowering of the electron-electron interaction 
through exchange by an amount /, here called the Stoner 
parameter. Because of the Pauli principle electrons with 
like spins are kept further apart and so their electrostatic 
Coulomb repulsion is, on average, weaker than for unlike 
spin electrons. This is the origin of Hund's rule as well 
as spin polarisation of itinerant electrons. Using these 
definitions of U and / we can write down E^ in terms of 



where the second derivatives are to be evaluated at the 
input density, i.e., m — 0. 

We also have, 

6p^ — — {Sp + Sm) , Sp'' — — [5p — 5m) 

from which we readily obtain the central result of this 
section. 



E^ = -USp"" - -ISp" - -Im\ 
2 2 4 4 



(5) 



Only the first two terms survive in the non spin polarised 
model described previously.^^'^^-^^ An associated expres- 
sion may be obtained from the LDA+J7 formalism as 
demonstrated in Appendix |^ Finally, we give the ex- 
pression for the tight-binding total energy including the 
magnetic terms, 



with 



and 



Etut = El + E2 



Ei=J2^^[p"Ho]+E^ 



(6) 



So 



R L L 



Ur - ^/r j Sq^ - ^/Rm^ 



in which p'^ is the spin density matrix, R labels atomic 
sites and Squ and V^{^ are as defined in equations (|B1[) 
and (|B2p in Appendix [B] There are no additional contri- 
butions to the interatomic force due to spin polarisationi^ 



III. RIGID BAND STONER-SLATER MODEL 



In the previous section we expanded the Hohenberg- 
Kohn total energy to second order around a non spin 
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polarised reference density. Alternatively one may ex- 
pand about a spin polarised density having a non zero 
trial magnetic moment.— We now show that in this case 
an expansion to first order is appropriate and that the 
resulting Harris-Foulkes functional may lead to the well 
known rigid band Stoner-Slater model,— i^i^li^ (usu- 
ally referred to as just the "Stoner model"). We re- 
call first that this is most readily illustrate d^°i'^^'^° us- 
ing the rectangular density of states, representing the d- 
band in a transition metal shown in figure [1] We imagine 
that majority spin electrons see an exchange and correla- 
tion potential lower than that seen by minority electrons 
by an amount proportional to the magnetic moment, 
m; the proportionality constant, /, being the "Stoner 
parameter." (Stoner uses the symbol a for this, / is 
Slater's usage.^) Then the rectangular bands are split 
by i^Ae = ±^Im and the change in band (kinetic) en- 
ergy due to magnetisation is 



■ Ae 



+ iAe 



FIG. 1: To illustrate the simple rectangular density of states 
model of ferromagnetism. A density of states which is con- 
stant and equal to g between band edges ±|VF is split 
by exchange into majority and minority spin densities by 
an amount proportional to the moment, m. (This can be 
achieved by flipping the spins of im electrons and realigning 
the Fermi levels to a common value.) By construction, we 
have Ae = m/g = Im. (After Pettifor,— figure 8.f2c) 



(pep rep pep \ 

/ gede + / gede — 2 / gede 

1 r2 2 
= --ff/ TO , 

using Ae = Im. In this estimate of the magnetic en- 
ergy the electron-electron interaction energy, —^Im? has 
been double counted, so it is subtracted to give 

Ai?.„, = J/to2 (1 - Ig) 

which is negative as long as Ig > 1, which is the sim- 
plest statement of the Stoner criterioni^ This particular 
model is pathological because AE,^^^ has no minimum as 
a function of to. This is a symptom of using a constant 
density of states, so that the kinetic energy is quadratic 
in to; that is, the fourth order term which is responsible 
for stabilising the ferromagnetic state is missing in the 
absence of structure in the density of states. 

According to Slater^ ferromagnetism arises from a com- 
petition between kinetic energy and on-site Coulomb 
electron-electron interactions. For an arbitrarily shaped 
density of states the kinetic energy increases compared 
to the spin-paired state when down-spin electrons are 
spin-flipped, since they must then be promoted into un- 
occupied states above the Fermi level. To develop a mag- 
netic moment, m, charge is transferred across the Fermi 
surface in small increments dTO, each increment costing 
more energy than the last as the down-spin states are de- 
pleted below the Fermi level and need to be taken from 
lower energy states and placed as up-spin electrons in 
higher energy states as these become successively occu- 
pied above the Fermi level. Generally speaking the larger 
the density of states near the Fermi level the smaller is 
the energy penalty involved. To counter this increase in 
kinetic energy there will be a decrease in energy due to 



a Hund's rule like exchange interaction and Slater— ar- 
gues that this takes the form —jIm?. Hence the total 
change in energy upon forming a magnetic moment to 

^2 . 12. 15. 32. 39 

Ai?„„, TO = - / — — - -Im\ 7 
2 Jo 5(™ ) 4 

which is clearly stationary at a generalised Stoner con- 
dition, namely Ig(m) = 1, where g{m) is the density of 
states averaged over the energy range spanned by flipping 
the ^TO spins; see figure 30 in the Varenna notes. ^''^ 

This is a rigid band model, requiring us to know only 
the non magnetic density of states. We can obtain an 
analogous expression for AE^-^^^ from a Harris-Foulkes 
functional, namely the first two lines of equation ([2|). In 
contrast to the second order theory in which the input 
density is non spin polarised, let us consider a trial den- 
sity which can be varied by changing its magnetic mo- 
ment while not affecting the total charge densityi^ We 
now have TOj,, = prj^ — prj and the trial Hamiltonian is 

where 

SO that a = +1 are the majority spins (ie, see a lower 
exchange and correlation potential). We now evaluate 
the first order total energy, 

i?W(TOj=^(V'^kl^'^ICk) 

- E / P^K^M - + E'Z + i?zz. 

When we compare this to its value when TOi„ = we 
obtain 

A£;_,(to„J = Ai;,,„,(TO.„) -f -^Iml (8) 
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FIG. 2: Contributions to the total energy in an LMTO''^ cal- 
culation for pure bcc-Fe relative to their values at m — 0. 
This is the Harris-Foulkes energy~2i^ -Ehf, as a function of 
the fixed magnetic moment of the input density, pi„ . Tc and r„ 
are core and valence kinetic energies and Eu is the Hartree en- 
ergy. We find that E^^ is almost exactly quadratic and hence 
its curvature is independent of m. Its curvature here is —0.04, 
giving / = 80 mRy compared to the value 65 mR y^'*''*^''^^ us- 
ing both Janak's method-- and that of Poulsen et al~ and 
68 mRy as calculated by Gunnarsson;^^ The fourth order term 
in Eu¥ which leads to a minimum at the observed moment 
comes from the kinetic energy. 



after evaluating the double counting in view of the fact 
that only the moment and not the density differ in the 
two cases, and using Ai?^" 
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As an illustration, we show in Figure [2] how a Harris- 
Foulkes energy varies with moment in pure bcc-Fe. Here, 
we have constructed an input density by superimposing 
free atoms^i^^ having a given magnetic moment so that 
the moment of the input density is a trial TOi,,. We then 
evaluate the Harris-Foulkes total energy functional and 
plot it against mi„. This is not exactly a rigid band cal- 
culation, but it serves to illustrate how the individual 
contributions to the energy vary with TOi„. In particular 
note that the kinetic energy increases, having both sec- 
ond and fourth order terms in mi„, while the exchange 
and correlation energy is found to be strictly quadratic. 
This is consistent with the Stoner picture and serves to 
show that the Stoner parameter / is independent of the 
moment and so may be taken as the same quantity in 
both equations ([6]) and ([8]) . Our estimate of / is of course 
not as good as a fully self consistent calculation as we in- 
dicate in the caption to figure [21 

We will use equation ([6]) to calculate density of states and 
total energy in sections |V] to IVIIII The rigid band pic- 
ture is particularly useful in interpreting complex mag- 
netic structures and arriving at an explanation of the 



enthalpy anomaly. Therefore in sections IIXI and |X] we 
employ equation ([H]) to find the total energy. 



IV. TIGHT-BINDING MODEL 



Our tight-binding model is specified by distance depen- 
dent matrix elements of the Hamiltonian and overlap, 
by Hubbard U and Stoner / parameters, and by a re- 
pulsive pair potential. We are motivated to employ the 
simplest possible scheme so as to maximise its predictive 
power relative to its complexityr^ Our starting point is 
the tight-binding theory of transition metals of Spanjaard 
and Desjonqueres*^ who propose a universal, orthogonal 
scheme in which Hamiltonian matrix elements have the 
form foe^'^'^ and the pair potential takes the form Be~^'^, 
where d is the bond length. These are intended to extend 
to nearest neighbours only in fee and hep metals and to 
second neighbours in the bcc structure. Spanjaard and 
Desjonqueres find a universal ratio p/q = 2.95 that fits 
well to the binding energy curve of Rose et al^ We have 
found this to be an excellent model for transition metals 
using an orthogonal basis of d-electrons^^ and adopting 
the canonical ratio for the three quantities /o, namely 

dda : ddir : ddS — — 6/o : 4/o : — l/o- 

Spanjaard and Desjonqueres provide values of the prod- 
uct qdo, where do is the equilibrium bond length, for 
most transition metals. Therefore the only adjustable 
parameters are /o which we adjust to the bandwidth 
calculated in the LDA, and the parameter B which we 
adjust to obtain the correct atomic volume (or lattice 
constant). This simple model having two adjustable pa- 
rameters then gives a good account of structural stability 
and elastic constants 

For a number of reasons, we wish to go beyond this very 
simple scheme in three respects, (z) We will extend the 
range of the exponentially decaying interactions; specifi- 
cally we encompass 58 neighbours in the bcc lattice. This 
has the attraction of employing an energy surface with- 
out discontinuities in a molecular dynamics simulation. 
Furthermore we have found this necessary to obtain a 
faithful reproduction of the LDA density of states, (ii) 
For this latter reason we also prefer to include s and p 
electrons in the basis, and (iii) to adopt a non orthogonal 
basis. We see a number of attractions from the inclusion 
of overlap which we discuss in Appendix [B] (see also the 
caption to figure 31 below). It is furthermore known that 
the neglect of sd-hybridisation leads to an overestimation 
of the magnetic moment of Fefiii^i^ Our procedure for 
obtaining the additional parameters is again motivated 
by simplicity and we adjusted the additional matrix el- 
ements to obtain a close comparison between the LDA 
and tight-binding density of states in bcc Fe. Thereafter 
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we merely adjusted fo to allow for the differences in d- 
bandwidth across the transition series. We use the same 
exponent in the overlap matrix elements as in the Hamil- 
tonian, but with a different prefactor, they thereby take 
the form soe~'^'^. We use g^Q^"^ = 3 for all dd interactions 
otherwise we set q — 0.5 bohr^^. We deviated from the 
canonical ratios in the non orthogonal case: 



dda : ddi: : ddS = — 6/o : 5/o : — 2.2/o, 

and furthermore used the ratio 

ppcr : ppir = 2 : — 1. 

We fix the on-site energy levels of the s and p atomic 
levels at 0.2 Ry and 0.45 Ry respectively, relative to the 
d-level. The remaining parameters are shown in table HI 




FeNM 



FeNM 



0.5 



Our values of the Stoner / are essentially those calcu- 
lated by Gunnarsson and others J^ii^i^'^i^ However we 
adjust these to obtain magnetic moments in agreement 
with the LSDA. 

Figure [3] illustrates the match between LDA and tight- 
binding densities of states in the non orthogonal and 
orthogonal d-only tight-binding models. Note that the 
canonical model is quite adequate in describing the es- 
sential features, namely the t2g {xy, yz, zx) bonding and 



{x^ - y 



r ) antibonding manifolds which sta- 



bilise the bcc structure at half band filling and the large 
density of states at the Fermi level, giep), which is re- 
sponsible for the ferromagnetic instability. To place the 
Fermi level exactly at the peak, it is necessary to choose 
the number of c?-electrons, N^, as an additional param- 
eter in the d-only tight-binding model; we set Nd = 6. 
However the three peak structure typical of bcc transi- 
tion metals and the smooth "U" -shaped pseudogap are 
less faithfully reproduced in the canonical model. 



FERRO- AND ANTIFERROMAGNETISM IN 
PURE IRON AND CHROMIUM 



It is quite clear that both canonical and spd tight-binding 
models predict ferromagnetism in Fe based in the Stoner 
criterion, Ig{eF) > 1, which in the simplest rectangular 
band models of Friedel^ and Pettifori2 is I /W > 1/5, 
where W is the width of the d-bandiiS In figure S] we 
show the self consistent tight-binding density of states 
compared to the LSDA. We find a self consistent mag- 
netic moment of 2.18/z^. The density of states of Cr 
is of course of practically the same shape as that of Fe 
but the Fermi level falls inside the pseudogap. In Pet- 
tifor's skewed rectangular c?-band theory ; ^'^1'^° antiferro- 
magnetism is predicted if 



FIG. 3: Density of states of non magnetic Fe using two tight- 
binding models. The upper panel shows the non orthogonal 
model having the parameters shown in table IH the dotted line 
shows the LDA density of states. The lower panel shows the 
density of states in the canonical d-band tight-binding model. 
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FIG. 4: Density of states in ferromagnetic Fe using the spd 
tight-binding model and compared to an LSDA calculation. 
The upper panel shows the minority and the lower panel 
the majority spins. Note that in an orthogonal tight-binding 
model, even using the fully self consistent scheme of sectionlTIl 
the two densities of states would be identical, only rigidly 
shifted. The inclusion of an overlap breaks this symmetry 
and it is seen here that this additional freedom acts signifi- 
cantly to improve the comparision with the LSDA. 



In this theory, the analogy is made between an AB binary 
alloy and an antiferromagnetic crystal having two sublat- 
tices, as does the bcc structure. In the alloy electrons will 
see a lower potential, say, at the A-site where the on-site 
energy level is lower than at the B-site by an amount Ae. 
In the common band model this leads to a skewing of the 
simple rectangular density of states, so that lower energy 
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TABLE L Parameters of our tight-binding model. Atomic Rydberg units are used throughout. 





ssa 


spa 


ppa 


sda 


pda 


dda 
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I 




fo So 


fo So 
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fo So 


fo So 






Cr 


-0.75 0.5 


0.5 -1.0 


1.0 -0.1 


-0.12 0.8 


-0.5 


0.18 




0.050 


Fe 


-0.75 0.5 


0.5 -1.0 


1.0 -0.1 


-0.12 0.8 


-0.5 


0.12 


340 


0.055 
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FIG. 5: Illustration of Pettifor's skewed rectangular band 
model of antiferromagnetism. The upper and lower figures 
show the densities of states on the two sublattices. A single 
one of these describes the situation in an AB alloy in which 
the electrons see a lower potential, say, at the A-atom whose 
density of states is accordingly skewed toward lower eigen- 
values as in the upper density of states of the top diagram. 
In the antiferromagnetic analogy, on each sub lattice the ma- 
jority spin electrons see a lower potential due to exchange 
interactions, they spend more time at that site and the den- 
sity of states is accordingly skewed. Spin up are the majority 
electrons at one sublattce, spin down at the other; hence the 
two diagrams, one for each sublattice. (After Pettifor,— fig- 
ure 8.12b) 
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FIG. 6; Density of states in antiferromagnetic Cr showing 
both the spd tight-binding model and an LSDA calculation. 
In the upper panels are shown both the up spins and the 
down spins from the lower panel reflected in the a;-axis to 
reveal the lower density of states compared to the up spin. 
The reverse situation pertains on the other sublattice. The 
lower density of states over the occupied density of states is 
the generaUsation of the skewed rectangular band picture. 



eigenvalues are generally associated with the A-sitc and 
vice versa. In this picture electrons in the lower energy 
single particle states spend more time at the A-site while 
overall charge neutrality is maintainedi^ In the antifer- 
romagnetic case (figure[S]) one says that up-spin electrons 
see a lower exchange potential at one sublattice and the 
down spin at the other. Each of their on-site energies are 
lowered through the exchange interaction (Hund's rule) 
by an amount As — Im, if / is sufficiently large, which 
favors aligned spins. Figure[B]shows that this effect is pre- 
dicted in the self consistent tight-binding model and com- 
pares the resulting density of states with the LSDA. The 
local antiferromagnetic moment m in the tight-binding 
model is predicted to be 0.74/1^^ in close agreement with 



the 0.70/jjj estimated from the LSDA spin density. 



VI. TRANSFERABILITY TO COBALT 



We begin discussion of energetics with the application 
of the Spanjaard and Desjonqueres model to Co. The 
approach we have taken is to adjust the parameter fo 
only to match the d-bandwidth of non magnetic bcc-Co 
calculated in the LDA. The resulting density of states is 
shown in figure [7] which also shows the density of states 
in hcp-Co to demonstrate the transferability of the band 
parameters to the observed structure of Co. 
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FIG. 8: Energy volume curves in Fe using the spd tight- 
binding model of table HI We used hep having ideal axial 
ratio. The dotted line shows non magnetic hep; the two curves 
merge as the antiferromagnetic moment vanishes with reduc- 
ing atomic volume. 



FIG. 7: Density of states in non magnetic Co: upper figure 
bcc, lower figure hep. Dotted lines show the LDA densities 
of states. Note that the parameters generated for the bcc 
structure transfer well to the observed hep structure. 



The remaining parameter, B, that enters the pair poten- 
tial was fitted to the calculated lattice constant of non 
magnetic bcc-Co. Table [II] shows the results of calcula- 
tions of both bcc and hep Co. The model is clearly re- 
markably predictive and argues strongly for the essential 
correctness of the Spanjaard and Desjonqucres approach 
coupled to the second order Stoner theory. Particularly, 
note that the tight-binding correctly predicts the stabil- 
ity of the hep over the bcc structure and also renders 
rather well the bulk moduli, both in magnetic and non 
magnetic forms. In connection with the Stoner / param- 
eter, we note firstly that the value, 68 mRy, quoted for 
the LSDA is not, of course, an input into the calcula- 
tion but this is the number calculated by other authors 
using the LSDA approach. Secondly, we note that we 
tried two values in the tight-binding model: / = 80 mRy 
gives a better value of the magnetic moment in bcc-Co, 
whereas this value gives a negative magnetic energy for 
hcp-Co thus predicting this phase to be non magnetic. 
Increasing / to 85 mRy corrects this but overstates the 
moment in bcc-Co. 



VII. PHASE STABILITY IN IRON 



We continue to look at the energetics by examining how 
the simple Spanjaard and Desjonqueres model describes 
the stability of the close packed structures in Fe. This 
has been addressed in detail recently;^ so for brevity we 
discuss only the bcc and hep structures at two atomic vol- 
umes, V/Vo = 1 and 0.88 where Vq = 11.82A3 is the ex- 
perimental atomic volume of bcc-Fe and the transition to 
hcp-Fe is observed^'^ to occur at about V/Vq — 0.88. Ta- 



ble [nT] shows that the predictions are less accurate than 
in the case of Co. We recall that very careful studies of 
the energetics in the LSDA have been made by Bagno et 
al.^"^ and by Stixrude et alM- The conclusions are that at 
V/Vq = 1; the most stable phase is ferromagnetic bcc-Fe, 
but that the energy volume curve for antiferromagnetic 
hcp-Fe intersects that for bcc-Fe and has a minimum at 
a lower energy at V^/Vb ~ 0.88. Hence the global predic- 
tion of the LSDA is that hep is the stable phase having 
a higher than ambient density. It is well known that this 
anomaly is removed by use of the so called generalised 
gradient approximation (GGA), although Bagno et al. 
point out that this is probably merely a coincidence aris- 
ing from the GGA favouring of both larger atomic vol- 
umes and larger magnetic moments as a general rule. As 
can be seen in figure [3 our tight-binding model rather 
closely follows the LSDA, but fails to reproduce the sta- 
bility of bcc-Fe even at the ambient atomic volume. Ta- 
ble mil shows also the predicted magnetic moments and 
bulk modulus. Note that we have used the ideal axial 
c/a ratio for hep at V/Vo = 1, but its measured value at 
V/Vfi = 0.88. 

Maybe it is not surprising that this very simple tight- 
binding model fails to describe the energetics of Fe. This 
is a very subtle problem even for the LSDA. The solution 
within tight-binding is rather simple however as has been 
demonstrated recently, and requires the use of a more 
complicated pair potentiali^ This is consistent with the 
observations of Bagno et al^ concerning the role of the 
GGA, and need not concern us further here, since in what 
follows we will discuss electronic structure and leave aside 
the question of structural energetics. 
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TABLE II: Energetic data for Co, eomparing tight-binding and LSDA ealeulations. Note that the only fitted values are the 
atomie volume of bee-Co, although the Stoner I has also been adjusted to agree with the LSDA moments. V/Vo is the atomie 
volume compared to experiment; Ai5h-b is the energy of the hep relative to the bcc phase; m is the magnetic moment; Ai5,„ag, 
the "magnetic energy" is the calculated energy difference between magnetic and non magnetic phases; K is the bulk modulus. 
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TABLE III: Energetics of Fe in the bee and hep crystal structures. Note that the tight-binding model incorrectly predicts that 
hep is stable at V/Vo — 1, but correctly reproduces the LSDA result that hep is stable at V/Vo ~ 0.88. The magnetic energies 
show that hep is only very weakly antiferromagnetic, especially at high pressure and this result is correctly reproduced by the 
tight-binding model. 
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VIII. ELECTRONIC STRUCTURE IN THE 
IRON-CHROMIUM ALLOY SYSTEM 



A. FeCr in the B2 crystal structure 



For the remainder of the paper we discuss the electronic 
structure of Fe-Cr alloys. It is very simple to construct 
a model for interactions between Fe and Cr by taking 
the geometric mean of the d-d hopping integrals and by 
moving the on-site d-orbital energies up and down by 
0.1 Ry. Thereby one would expect a small charge transfer 
from Cr to Fe, since the latter is more electronegative. 
To control this charge transfer we apply a Hubbard U of 
1 Ry. Our model deviates in this way slightly from the 
usual ansatz of local charge neutralityi^ 

The B2 alloy FeCr has a positive heat of formation and 
hence does not exist i^^i^^ Nonetheless it presents an inter- 
esting case in which to discuss the competition between 
ferro- and antiferromagnetism. One might expect this 
ordered alloy to be antiferromagnetic since the Cr sub- 
lattice could prefer to align antiferromagnetically with 



the neighbouring Fe atoms. But the non magnetic den- 
sity of states clearly shows a large density of states at 
the Fermi level and one expects the Stoner criterion to 
apply and lead to ferromagnetism. However, it turns out 
in the tight-binding model that both ferro- and antiferro- 
magnetic solutions can be found depending on the value 
of the Hubbard U ; but in the physically correct limit of 
large U the alloy is ferromagnetic in agreement with the 
LSDA. To begin with, figure [9] shows Fe and Cr atom 
projected densities of states in non magnetic FeCr. We 
observe that the small amount of charge transfer permit- 
ted by the self consistent tight-binding leads to a closer 
agreement with the LDA than the non self consistent 
tight-binding density of states. 

Figure [TO] shows the density of states in the self consis- 
tent spin polarised tight-binding calculation employing 
a Hubbard t/ of 1 Ry. The result is in close agreement 
with the LSDA. The local moments on the Fe and Cr 
are 1.14 /l/^ and 0.71 /^^, in reasonable accord with the 
estimated local moments in the LSDA, namely 1.46 ^ 
and 0.34 Figure [TT] shows the local moments as a 
function of the Hubbard U where we find an unphysical 
regime if charge transfer is allowed to occur. In that case 
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FIG. 9: Atom projected densities of states in non magnetic 
FeCr. The top panel shows a non self consistent tight-binding 
calculation in which the charge transfer is seen to be some- 
what smaller than that found from the self consistent tight- 
binding calculation in the centre panel. This self consistent 
result agrees better with the LDA in the lower panel. 



we find an equal number of electrons in the spin up chan- 
nel, while in the spin down there is a larger population 
on the Fe than the Cr site leading to antiferromagnetism. 



B. Chromium as a dilute impurity in iron 



Whilst FeCr is ferromagnetic, a Cr atom in dilute con- 
centration in Fe becomes antiferromagnetically ordered 
with respect to the Fe host atoms>2^ We find that the self 
consistent tight-binding model reproduces the LSDA re- 
markably well in detail, and furthermore offers an expla- 
nation rather more readily than the LSDA. We illustrate 
this using a unit cell of 16 sites in the bcc Fe lattice, in one 
site of which an Fe atom is replaced with a Cr atom. In 
figure [T21 we show local densities of states projected onto 
the Cr and its neighbouring Fe atoms, both using LSDA 
and tight-binding. Note how the local density of states 
projected onto the Fe atoms neighbouring the Cr impu- 
rity is hardly different from that of bulk Fe. It is curious 
that the Fe does not accommodate itself to the presence 
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FIG. 10; Atom projected densities of states in ferromagnetic 
FeCr. The self consistent tight-binding model is in close agree- 
ment with the LSDA. 
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FIG. 11: Local magnetic moments (in /i^) on the Fe and Cr 
atoms of B2 FeCr in the self consistent tight-binding model as 
a function of the chosen Hubbard U parameter. For compar- 
ion the estimated local moments from the LSDA are shown 
as horizontal lines (of course, the Hubbard U is "built-in" to 
the LSDA and cannot be varied. This may be regarded as an 
advantage and a disadvantage: in the tight-binding one may 
observe the role of parameters such as / and U by varying 
them) . Note the transition from antiferro- to ferromagnetism 
as U is increased. The limit of local charge neutrality leads 
to rather smaller moments than found in the LSDA, but does 
predict the correct magnetic ordering. 
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FIG. 12: Densities of states in an ordered FeisCr alloy on the 
bcc lattice, (a) and (c) show the local density of states pro- 
jected onto the Cr atom, respectively using LSDA and tight- 
binding. The minority spin density of states in the lower panel 
is repeated, using a dotted line, for comparison by reflection 
about the abscissa, (b) and (d) show the local density of states 
projected onto the Fe atoms neighbouring the Cr, again using 
LSDA and tight-binding respectively. A dotted line shows the 
density of states in bulk Fe. 



of the Cr impurity. On the other hand, the Cr projected 
density of states is greatly perturbed from its bulk, as 
may be seen by comparison with figure [SI The most 
prominent feature is a narrow resonance in the occupied 
majority spins, which is almost completely unhybridised 
with the neighbouring Fe minority spins. We show in 
figures [T3] and [T3] the densities of states from figure [T^] 
projected into the t2g and Cg manifolds. It becomes clear 
that this prominent feature arises from strongly localised 
states of xy, yz and zx character. 



IX. THE MAGNETIC RIGID BAND MODEL 



It is clear from a comparison of figures [T^] and [S] that a 
rigid band approximation would be a very poor descrip- 
tion of alloying in the Fe-Cr system. A recent calculation 
using the coherent potential approximation in the LSDA 
has been made,-- but in this case, the densities of states 
do not very well resemble those shown here in figure 1121 



FIG. 13: LSDA local densities of states in FeisCr projected 
into the t2g and eg manifolds. This shows the non bonding 
resonance on the Cr impurity to originate from the xy, yz, zx 
and d-orbitals. 



However, we may use the rigid band model described in 
section Hm in which the input density is constructed hav- 
ing a trial moment. Indeed as seen in figure 1141 such 
a trial density (to be described in detail below) gives a 
very faithful reproduction of the self consistent density of 
states. In the simplest example, that of the non magnetic 
density of states of Fe shown in figure [31 a plot of AE,^^^ 
from equation ([8]) versus m is shown in figure I15i hav- 
ing the characteristic double- well^ structure with minima 
at m = and a magnetic energy of 21 mRy; these 

values may be compared with those from the self con- 
sistent tight-binding calculation in table [TTTl viz. 2.18/ig 
and 17 mRy. The small discrepancies arise from the self 
consistent calculation allowing the shape of the spin den- 
sities of states to be different from the input, non self 
consistent densities. As mentioned in the the caption to 
figure [31 above, this is entirely due to the use of a non 
orthogonal tight-binding model. 

We can now use this simple construction to interpret the 
stability of the antiferromagnetic alignment of the Cr im- 
purity in Fe. A trial spin polarised density is constructed 
by imposing a moment of -f 2.2/^ on each of the Fe atoms 
and a trial moment m on the Cr impurity. The associated 
bandstructure energy difference is found to which jlrri^ 
is added, in which we take / = 50 mRy from table[Tl The 
magnetic energy plotted against m is shown in figure 1161 
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FIG. 16: Magnetic energy versus trial local Cr impurity mo- 
ment in FeisCr. Note there is only one, antiferromagnetic 
solution. 
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FIG. 14: As figure 1131 but using the tight-binding model. 
In the upper panel a dotted line shows the density of states 
arising from a trial input non self consistent consistent density 
with imposed moments on the Fe and Cr atoms (see the text). 
Note, firstly, the excellent detailed agreement with the LSDA 
in figure [12] and secondly the close similarity between the self 
consistent and non self consistent densities of states. 
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FIG. 15: Magnetic energy versus trial moment in the rigid 
band Stoner model for bcc-Fe (compare with figure[2|). Curves 
of this type were first computed by Slater. 



Only one, antiferromagnetic, solution is found, having a 
local Cr moment of 2.42/^^, which is close to the moment 
of 2.37 fig found in the self consistent tight-binding calcu- 
lation. (The estimated Cr local moment from our LSDA 
calculation is 2.08fig.) Although there is no ferromag- 
netic solution, it is instructive to plot the trial densities of 
states for trial local moments of 2.37fig in both antiferro- 
and ferromagnetic alignments. These are shown in fig- 
ure [iTl Neither looks at all like the density of states of 
pure Cr in figure [51 this is because to develop an anti- 
ferromagnetic state requires the cooperation of two sub 
lattices, which cannot be achieved by isolated Cr atoms 
or small clusters of these (say, fewer than nine atoms) in 
a bcc-Fe host. This is why it is the Cr density of states 
that has to accommodate itself to the underlying Fe elec- 
tronic structure, and this lies at the heart of understand- 
ing the enthalpy of mixing and the phase diagram in the 
Fe-Cr system. Figure [17] helps to explain why isolated 
Cr impurities do not align themselves ferromagnetically 
with the host Fe. To do so would require a density of 
states essentially that of pure ferromagnetic Cr, and this 
phase is unstable with respect to the observed antiferro- 
magnetic phase in Cr. The alternative is to align antifer- 
romagnetically, and this causes the density of states to 
adopt a shape quite unlike that in pure Cr while at the 
same time there is a complete lack of cooperation from 
the very stable bcc-Fe density of states, practically the 
same as pure Fe even on the Fe atoms neighbouring the 
impurity. 
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FIG. 17: Trial densities of states constructed by imposing 
fixed magnetic moments onto the Cr impurity in FeisCr. As 
already seen in figure 1141 the trial density for the antifer- 
romagnetic alignment is very close to the corresponding self 
consistent density of states. Because there is no energy mini- 
mum at ferromagnetic alignment (see figure lT^ such a density 
cannot be achieved in a self consistent procedure. This illus- 
trates the usefulness of this construction. 



X. ORIGIN OF THE REPULSION BETWEEN 
CHROMIUM IMPURITES AND THE 
ENTHALPY ANOMALY 



Now we ask what is the stable magnetic structure of two 
Cr impurities placed as nearest neighbours in Fe? We 
go straight to the predictions of the tight-binding model 
shown in figure [THl We make trial spin densities having 
the Cr spins parallel or antiparallel to each other and plot 
the magnetic energy as a function of their moment. In 
the case that they are antiparallel and assuming the two 
moments to have same magnitude, we find a double well 
as expected. The more stable structure is for both spins 
to be aligned parallel to each other, but to be antiferro- 
magnetically aligned with the spins of the Fe host. In 
fact the antiparallel state is unstable and we find that if 
the constraint is removed in a self consistent calculation 
this reverts to the parallel state. 

If the Cr atoms are placed at second neighbour positions, 
with their spins aligned parallel to each other, we find an 
energy versus magnetic moment very similar to that of 
the single impurity in figure 1161 In fact our LSDA and 
tight-binding calculations (not presented here) show the 
densities of states and magnetic moments to be very sim- 
ilar in these two cases; indeed the LSDA local Cr moment 
is a little larger in Fei4Cr2 than in FeisCr as seen also in 
the tight-binding model by comparing figures and 1181 
This latter figure now illustrates rather clearly the origin 



of the repulsion between Cr impurites in bcc-Fe. The 
energy is lower when the atoms are placed at next near- 
est neighbour positions as long as spin polarisation is 
allowed; otherwise the energy ordering is reversed as is 
also found using LSDA calculations.^^ Furthermore since 
the jIm? term is the same in both cases this is clearly a 
handstructure effect. 

We can now offer a more detailed explanation for the 
anomaly in the enthalpy of mixing of Cr in Fe. In most 
of the concentration range Cr prefers to cluster together 
to allow sufficient atoms to cooperate towards providing 
the two sublattices required to establish the antiferro- 
magnetic state. Hence the enthalpy of mixing is positive 
and spinodal decomposition is observedi^ Conversely at 
low concentrations, the Cr may appear as isolated impu- 
rities stabilised by the change in spin polarised density of 
states which has quite a large weight at the bottom of the 
band as seen in figure 1171 These isolated impurites repel 
each other, as already found by Klaver et al.^ shown 
clearly in our figure [TSl hence at low concentrations the 
enthalpy of mixing is negative, but only while the concen- 
tration of Cr is sufficiently low for the Cr-Cr repulsion 
to dominate. Our present modelling explains the nearest 
neighbour repulsion in detail. The LSDA calculations^^ 
also showed that the Cr-Cr repulsion extends to second 
neighbours and beyond, these longer ranged interactions 
contribute significantly to the total repulsive energy of 
a pair; furthermore they are present even when the sys- 
tem is forced to be non spin polarised, when the near- 
est neighbour repulsion collapses. An explanation of the 
longer ranged repulsion remains to be found in the band- 
structure. 



XL CONCLUSIONS 



We have described how itinerant magnetism can be in- 
corporated into our self consistent polarisable ion tight- 
binding model. This results in an additional parame- 
ter, the Stoner /, which we identify as minus twice the 
curvature of the exchange and correlation energy as a 
function of magnetic moment. A first order expansion of 
the Hohenberg-Kohn functional leads to the rigid band 
Stoner-Slater model. We show that a very simple pa- 
rameterisation of the tight-binding model is possible that 
gives a faithful reproduction of the energetics and elec- 
tronic structure of the LSDA. The parameters of the 
model are easily transferable between the first row tran- 
siton metals and their alloys. The simplest form of pair 
potential is quite adequate, except in the case of Fe, to 
reproduce structural stability and bulk modulus. Armed 
with this model we address outstanding questions re- 
lated to solution and clustering of Cr impurities in bcc- 
Fe. The rigid band, fixed moment approach proves to 
be very useful in reproducing LSDA results and predict- 



14 




FIG. 18: Magnetic energy versus trial local Cr impurity mo- 
ment in Fei4Cr2, having the two Cr atoms as nearest neigh- 
bours, NN, and next nearest neighbours, NNN. The solid line 
denotes the energy of the pair of Cr atoms having their spins 
aligned parallel to each other. The broken line refers to the 
two Cr spins being aligned antiparallel to each other. The 
dotted line is the energy in the case that the two Cr atoms 
are separated as next nearest neighbours. The zero of energy 
in this plot is the energy of the non magnetic NN case. Hence 
the graph shows the lowering of energy as a result of mov- 
ing the Cr atoms from NN to NNN positions as long as spin 
polarisation is permitted. The repulsion between Cr impuri- 
ties is thereby revealed firstly as a bandstructure effect and 
secondly as an effect of the magnetism. 



ing magnetic structure and energy of complex transition 
metal alloy systems. This provides a powerful framework 
within which to explore complex magnetic structures in 
transition metals generally. The model is based in the 
correct physical picture, namely itinerant magnetism re- 
sulting from a competition between kinetic, or band, en- 
ergy described by inter-site one electron hopping matrix 
elements of the non self consistent tight-binding Hamilto- 
nian; and on-site exchange and correlation parameterised 
through a single Stoner parameter. Because the tight- 
binding approximation is particularly simple and trans- 
parent we believe that this approach will find a number 
of applications in this area in the future. 
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APPENDIX A: CONNECTION TO LDA+C/ 

We may arrive at an expression similar to (O from the 
starting point of the theory of LDA+Ui^ The usual no- 
tation is to write n'^ for the number of electrons or occu- 
pation number in, say, a d-band with quantum number 
m (not to be confused with the magnetic moment) and 
spin a. Then defining U and J as spheridised, orbital in- 
dependent Coulomb and exchange integrals, the on-site 
electron-electron interaction energy igS^.sSjSQjeOjei 

E'' = lu nl, n-r + \{U^J) J] , 

mm' (7 Tm-n' a- 

= \up' - \jY.{n' - \{u ~ j)Y^{Kf . 

(7 ma 

(Al) 

The first line shows in its first term unlike spins inter- 
acting through the Hubbard U , and in the second term 
like spin electrons interacting through a Hubbard term 
reduced by an amount J as explained at the end of sec- 
tion [Hi This term explicitly requires m ^ m' in the sum: 
as two electrons cannot occupy the same state accord- 
ing to the Pauli principle this would otherwise give an 
interaction between an electron and itself. Hence the on- 
site electron-electron interaction properly includes the so 
called self interaction correction present in Hartree-Fock 
theory but not in the LSDA»22 The second line^S. follows 
directly after some algebra, expressing 



The three terms resulting in the second line of (|Al|l are 
respectively a direct Coulomb term, an exchange term 
and a term which is of lower order of magnitude com- 
pared to the first two and which would amount to admit- 
ting an orbital dependent potential. In the spirit of the 
LSDA we neglect^'^ this last term and by differentiation 
we find for the potential seen by an electron with spin a 
as a result of electron-electron interaction. 



Up - Jp° 



and so the exchange splitting between up and down 
spin energy levels is approximately Ae ^ — = 
— J{p'^ — p~) = —Jm. After some further algebra again 
neglecting the third term in (|A1|) we may also write 



E^ = Uu-\j)p'-\jm' (A2) 



which is equivalent to our expression ([5]) for E2 in sec- 
tion [ill after identifying the exchange integral J with the 
Stoner parameter /. Note, however, that E^ is not an 
energy to second order in any charge density difference. 
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but it could be cast into such a form if we make an ex- 
pansion of the total energy in a generalised mean field 
multiband Hubbard model. We wish to emphasise two 
points here, (i) Both exchange and correlation are con- 
tained in equations ([5]) and (|A2p . the effective Coulomb 
integral being reduced to U—^J by exchange. Indeed it is 
well known that the exchange-only Kohn-Sham-Gaspar 
potential gives a poor description of itinerant magnetism 
by overestimating the tendency to magnetism in transi- 
tion metals. (ii) As in LSDA, equations (O and (K2\i 
are functionals of the spin density only and lead to or- 
bital independent potentials. It is clear, though, from 
the foregoing how to recover the self interaction correc- 
tion (at least in on-site terms in the Hamiltonian) in a 
tight-binding context in which the potential seen by an 
electron is orbital dependent. 



APPENDIX B: NON ORTHOGONAL SELF 
CONSISTENT TIGHT-BINDING 



There is a number of benefits of adopting a non orthogo- 
nal tight-binding basis. It is widely believed to result in 
a more transferable model. In addition it admits the con- 
cept of bond charged As we now demonstrate this allows 
the self consistency to adjust the hopping integrals as well 
as on-site matrix elements of the Hamiltonian. We re- 
call that our self consistent polarisable ion tight-binding 
modeiiSiii^'^ is couched in terms of multipole moments 
of charge with respect to neutral, spherical atoms having 
valence electrons. Hence the self consistent charge 
transfer to a site labelled by its position R in units of 
electron charge, e, is 

<5(7R = - = Qro. (B1) 

Higher moments of the charge develop as a result of crys- 
tal field splitting and these are denoted Qrl, in which L 
is a composite index subsuming both angular momenta: 
L = {£m}. The Madelung potential (energy) at site R 
due to multipoles at sites R' is 

R'/R L' 

B is a generalised Madelung matrix )^^'^^ related to the 
structure constants of LMTO theoryji^ For monopolc in- 
teractions, we write 

Boo (R'-R) = p^^^C/rr.. 

The transfer of charge is resisted by a "Hubbard poten- 
tial," 

- ?7rQro. 

In the orthogonal self consistent tight-binding model, 
these potentials are used to adjust the on-site matrix 



elements of the Hamiltonian, both on-site energies and 
off-diagonal crystal field terms. The increments to the 
Hamiltonian are 

VrlrL' = V-^SlL' + E ^i"^«'^"C'LL'L" 
L" 

in which Cll'L" are the Gaunt coefficents that enforce 
the selection rules and Auiin are new parameters con- 
trolling the strength of the crystal field splittingj^iiiSii^ 
These may be adjusted, for example, to reproduce crys- 
tal field splittings in ab initio bandstructures or dipole 
moments in molecules. 

If we include an overlap matrix S'rlr'L', then solving 
the generalised eigenproblem leads to normalised eigen- 
vectors and the charge at site R is 

gR = ^E/"k E ((^rI^rlr'L'C^II + c.c.) 

nk R'L'L 

Ef \ " l/^nk 1 2 
Jrik \ ^'R.L\ 

nk L'L 

+ 2^^In\!i E (^RL^R.iR'i'C'RL + C-C.) . 
nk R'L'L 

Here, a bar and "c.c." imply complex conjugation. /„k 
are occupation numbersSS of the state at wavevector k 
and band index n, as used say in Fermi-Dirac or gener- 
alised Gaussian Brillouin zone integration;^ or the lin- 
ear tetrahedron method.— The final term amounts to a 
bond charge which is absent in orthogonal tight-binding 
models. To extract the bond charge explicitly, we have 
defined O = S — 1 and since the norm is conserved sepa- 
rately at each k-point, we work with Bloch transformed 
matrices, such that, for example, 

'S'rlr'L' = E '5'r+tlr'L' e* "^, 

T 

where T are the translation vectors of the crystal lattice. 

For simplicity we allow the overlap to make contributions 
only to the monopole moments of the charge; higher mo- 
ments are defined as in the orthogonal case so that for 
£ > we have^iiii^iiaiSS 

QrL = E-^"k E ^RL'C'rl"^«'^"C'lL'L"- 
nk L'L" 

We now find increments to the hopping integrals as a 
result of the self consistent redistribution of bond charge. 
These areii 

^RLR'L' = ^Pr + ^R')O^LR'L' 

where 

R'O- 

R' 
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FIG. 19: Hubbard and Madelung contributions to the force 
on atom R. Circles are intended to represent changes in 
monopoles arising from the displacement of the atom at R, 
which modifies the charge on both sites R and R' through 
the scaling of the overlap matrix elements with bond length. 

To preserve the norm these need to be updated directly 
into the Bloch transformed Hamiltonian. Note that only 
monopole terms enter here as a result of our definition 
of the higher multipoles without reference to the overlap. 
Dr is the sum of Hubbard and point-charge Madelung 
potentials at site R. 

There are also new terms in the interatomic forces. Ac- 
cording to the Hellmann-Feynman theorem the force is 
obtained from the derivative of the energy, taken while 
keeping the wavefunction frozen. In an orthogonal tight- 
binding model multipole moments do not change under 
this constraint when the atom at R is displaced; hence 
the only contribution to the force from self consistent, 
second order terms in the energy is the classical electro- 
static term, 

■pes _ 1 2 ^ OBl'L" (R" - R) ^ 
^R--2 Qr'L' Qr"L"- 

R'^R" 
L'L" 



However in a non orthogonal model, even at fixed eigen- 
vectors, displacement of an atom will lead to changes in 
the bond charges with its neighbouring atoms as a result 
of the changes in the overlap matrix elements. There are 
two new contributions to the interatomic force. Since we 
are concerned with derivates of the overlap matrix, we 
will require the quantity 

,RLR'L' 

LL' 

where 

^Pr.RLR'L' = ^^Pr',RLR'L' 

1 r ( /^nk ^'S'rlR'L' /^Jik I ^ ^ ^ 

Then for the first contribution we find 

FK = -E(^K+^R')5P^R' 

R' 

and for the Madelung contribution, 

R' 

is the Hubbard potential, and V^q is the £ — 
component of the electrostatic potential (jB2p seen at R. 
These two contributions to the interatomic force are open 
to a quite simple interpretation if we make reference to 
figure [m 

When the atom at R moves, its own monopole moment 
changes by virtue of overlap with an atom at R'. This 
leads to a change in Hubbard potential (energy) at site R 
and hence a force (figurefTWa')'). This change in monopole 
moment at R will result in a modified electrostatic inter- 
action with a multipole moment at a third site R" (in- 
cluding the possibility R" = R') described by the matrix 
element i?oL" (R" ^ R)- This leads to the first Madelung 
contribution, shown in figure [TW c). The same movement 
also induces a change in the monopole moment at site R' 
giving rise to the second Hubbard contribution, shown in 
figure flW b). The second Madelung contribution, illus- 
trated in figure [TWd). corresponds to the force associ- 
ated with the electrostatic interaction between a multi- 
pole at R" (admitting the possibility that R" = R) and 
the modified charge at R' through the Madelung matrix 
element Bql" (R" - R')- 
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